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^ , model for vortex line fluctuations in a high Tc superconductor in an external magnetic field. A density of 
vortex lines of / = 1/25 is considered. We find two sharp phase transitions. The low T superconducting 
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entangling, cutting, and loop excitations. An intermediate phase is found, which is characterized as 
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superconducting properties is the direction parallel to the magnetic field, but normal behavior in planes 
perpendicular to the field. A detailed analysis of the vortex structure function is carried out. 
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I. Introduction. 

While an understanding of the microscopic mechanism leading to high temperature super- 
conductivity has remained elusive, much recent research has focused on obtaining a phenomeno- 
logical understanding of the behavior of these materials in the presence of applied currents and 
magnetic fields. Such treatments have been based on the familiar Landau- Ginzburg approach 
for type II superconductors, which assumes only that macroscopically, the superconductivity is 
adequately described by a complex order parameter, as is the case for classical BCS materials. 
In this context, the high materials are generally believed to differ from the classical supercon- 
ductors, in that fluctuation effects are dramatically enhancedfl], |^, ||. The relatively high values 
of Tc, anisotropy, and ratio of magnetic penetration to coherence lengths, k = X/C,o, in these 
materials lead to a significantly larger critical region in temperature, about Tc. 

In the presence of an applied magnetic field Hd < H < Hc2, the low temperature phase of a 
pure system is the Abrikosov vortex line lattice. As temperature increases, it has been argued 
that fluctuations will cause this vortex line lattice to melt at a temperature well below that given 
by the mean field Hc2(T) line, leading to a vortex line liquid phase |, ||, |^. Experimental 
evidence for such a scenario has been cited in mechanical oscillator measurements 0, as well 
as by the observation that the onset of sizable reversible magnetization generally occurs at a 
temperature well above that where resistivity vanishes[0], an indication of strong fluctuation 
effects. Several theoretical works have sought to describe the nature of this vortex line lattice 
melting transition, as well as predict properties of the vortex line liquid phase. One commonly 
used approach is an effective elastic theory [Q, H, ^, in which the energy of transverse fluctu- 
ations of the vortex lines from their lattice positions, is computed in terms of non-local elastic 
constants determined from the mean field Landau-Ginzburg theory [Q. The Lindemann criterion 
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is then used to estimate melting, when such transverse fluctuations become a sizable fraction 
of the vortex line spacing. Another theoretical model has been the "two dimensional boson" 
approximation 0, 0, in which vortex lines, fluctuating in a directed fashion along the direction 
of the applied magnetic field, are viewed as the "world lines" of two dimensional interacting 
bosons. Many-body methods applied to the two dimensional boson problem have then been 
used to estimate melting in the dilute line limit H, as well as to compute properties of the 
vortex line liquid 0. 

The above theoretical models suffer from several important limitations. The effective elastic 
theory loses validity once the vortex lines are in a liquid phase. The "2(i boson" approximation 
assumes a simplified vortex line interaction, that applies mainly when lines do not wander too 
much from parallel alignment with the applied magnetic field. Both models assume that the 
vortex lines thread the system in a directed fashion; each line passes through each plane perpen- 
dicular to B only once. Closed vortex line loop excitations, which have been argued to play a 
significant role in the B = case]^, |I^], as well as in the finite B case|]rT|, are ignored. In view of 



this, it is useful to have a well defined statistical model, for which numerical simulations can be 
carried out, and in which the effects of vortex line fluctuations may be quantitatively examined. 

In this paper we present the results of such a study, based on a simplified model which 
we believe captures the essential physics of a fluctuating Landau-Ginzburg superconductor; the 
uniformly frustrated three dimensional XY model. We have previously introduced this model in 



two earlier works |T2|, [T3|] . Here we provide details, as well as surprising new results. We find that 
as the vortex line lattice is heated, the system undergoes two distinct thermodynamic transitions. 
First there is a transition characterized by the loss of phase coherence in planes perpendicular 
to the applied magnetic field, while phase coherence is retained parallel to the field. In this 



intermediate region, the vortex line lattice disorders in the transverse direction, while individual 
lines remain essentially straight and unentangled from each other. Then at a higher temperature, 
phase coherence is lost in the direction parallel to B as well. This high temperature region is 
characterized by vortex line wandering and many closed loop excitations. The vortex lines are 
highly entangled, yet they easily cut through one another. 

The remainder of our paper is organized as follows. In section II we present a detailed 
description of our model, and the approximations involved. In section III we present the results 
of our Monte Carlo simulations: section IIIA discusses computation of the helicity moduli, which 
measure the loss of phase coherence; section IIIB discusses the nature of vortex line fluctuations 
in the various thermodynamic phases; and section IIIC presents a detailed analysis of the vortex 
structure function, which is used to extract effective elastic constants and correlation lengths. 
In section IV we summarize our conclusions. 
II. Description of the Model 

It is generally believed that the phenomenological behavior of the high temperature super- 



conductors is well described by the familiar Landau- Ginzburg free energy functional [|14|, 

Fim, m] = J d'r + i/^lV^r + + + i-|V X ipj (1) 

— * 

where ip is the superconducting order parameter, A the magnetic vector potential, and $o = 
hc/2e the flux quantum. The local magnetic induction is 6(r*) = V x A{f). The first two 
terms represent the condensation energy, the third term is the kinetic energy of the fiowing 
super currents, and the fourth term is the magnetic field energy. We choose to work throughout 
this paper with the Helmholtz free energy, which is a function of the average magnetic induction, 
B = (l/Vol) J d^rb{r), inside the bulk superconductor. The applied magnetic field H, may 
in principle be determined by an appropriate thermodynamic derivative of the Helmholtz free 



energy, but its value is of no direct concern for our calculations. For simplicity, we consider here 
only a material with isotropic couplings, and indicate later how the effects of anisotropy may be 
included in our model. 

The mean field solution, appropriate for low temperature superconductors, is obtained by 
minimizing Eq.(|l|) with respect to both ip and A, subject to the constraint of constant B. For 



B = 0, one gets the uniform solution ^jj{f) = y\a\/P = ipo in the superconducting phase, a < 
{ipQ = in the normal phase, a > 0). For < -B < if, in the mixed phase of a type II 



superconductor, the result is the familiar Abrikosov flux tube lattice [|14[. Each flux tube carries 
one $0 of flux, and consists of a vortex line in the phase 9(r) of ip{r) = \^{r)\e^^^^ . The vortex 
line is surrounded by a normal region of radius the coherence length. 



eo = h/^2m\a\ (2) 
and the magnetic field is confined to within a region of radius the magnetic penetration length. 



A = ^JmcyWne^ltpol'^ (3) 
The density of flux tubes is B/^o, hence the separation between vortex lines is 



%/B (4) 

For the high temperature superconductors, we now wish to compute thermodynamic quantities, 
averaging over all fluctuations of ip and A (subject to the constraint of constant B), weight- 
ing conflgurations with the Boltzman factor q-^'I'^^^I/'^bT _ order to carry out this average 
numerically, with Monte Carlo simulations, we make the following approximations. 

For A ^ a^, the magnetic fields of the individual fiux tubes strongly overlap, leading to an 
approximately uniform local magnetic induction. In this limit, we therefore make the approx- 
imation of replacing the fiuctuating A with a fixed one, giving a spatially uniform magnetic 



induction B = V x A. We will discuss the consequences of this approximation on the effective 
vortex line interaction in more detail below. Next we discretize[^ our model by placing it on a 
periodic cubic mesh of lattice constant a, resulting in the free energy functional 

F[^,,A,] = a^j: |(a + ^)|^/ + Im'] - i^E^K^^'^-^^ + cc. (5) 
j I ) (ij) 

where we have dropped the now constant magnetic energy term. Here i, j denote sites of the 

cubic mesh, (ij) denote nearest neighbor bonds of the mesh, and 

A,^ — I A-dl (6) 



$0 

Finally, we choose the lattice constant of the numerical mesh a = the coherence length of 
Eq.(0). Thus the decay to zero of the superconducting wavefunction ip at the normal core of a 
vortex line, takes place within one unit cell of the mesh. Outside the normal core, we can make 
a London approximation and set = ipoy a constant. Thus on the mesh, we consider only 
fluctuations of the form ipj = ip^e^^^ , where the amplitude is fixed, and only the phases 6j vary. 
Details concerning the normal cores of vortex lines are treated approximately by the small length 
cutoff a = ^0 imposed by the mesh. Applying this to the free energy in Eq.(^ results, within 
additive constants, in the effective Hamiltonian for a frustrated three dimensional XY model. 



F[e,] = -JoJ2cos{e,-e^-A,) (7) 

(ij) 



where Eqs.(0,|5D, $o = hc/2e, and = a give. 



- 16^ 



Anisotropy, or impurities, could now easily be incorporated into the model, by letting the coupling 
Jo vary appropriately on different bonds of the mesh. In this work we continue to consider only 
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the case of isotropic couplings. We expect that the qualitative behaviors we find, will be to a 



large extent the same as in anisotropic systems [p!2|, |T6| 



The Hamiltonian above has been studied extensively in two dimensions as a model for 
Josephson junction arrays[0, Il8|, and in two and three dimensions as a model for granular 



superconductors [0]. Its application to behavior in the high temperature superconductors, is best 
realized in the limit ^ c^i; ^ ^i where the above approximations should apply. This condition 



may be rewritten as Hd <^ H <^ Hc2- For the high Tc materials, where k = A/^o ~ 100 is 
large, this will cover a large range of experimentally accessible magnetic fields H. For YBaCuO, 
for example, this condition holds for fields on the order of a few Tesla. 

Much of the theoretical work concerning fiuctuations in the high Tc superconductors, has 
been discussed in terms of fiuctuations of the vortex lines induced by a finite B. It is therefore 
of interest to consider the effective vortex line interaction in our model (0). Towards this end, 
it is convenient to replace the cosine function of Eq.(0) with the periodic Gaussian, or Villain 
function [gg. 

Fm = JoY.v{e.-e,-A,,) (9) 



where 



oo 1 

Via) = -(T/Jo) ln{ ^ exp[-- Jo(« - 2nmf /T]} (10) 



m=— oo 



This has the effect of eliminating the coupling between spin wave and vortex excitations of the 
phase variables 6i in the Hamiltonian (|^). 

Making a standard duality transformation [21] on this new Hamiltonian (j^), gives the effective 
Hamiltonian in terms of the vortex degrees of freedom. 



F^nM = 27rVoEEK(i) " fMM) ' fMG{n-r^ (11) 

ij M 



Here is the integer vorticity in the phase 9, going around the plaquette with normal fi = x, 
y, z of the cubic unit cell centered at the dual mesh site i, and /^(i) is the number of flux quanta 
of magnetic induction through this plaquette, 

2^^(i) = '^f^-d^= ^Ok + + Aim + Arnj (12) 

where the sum is taken counterclockwise over the bonds forming the edges of the plaquette. For 
our uniform magnetic field approximation, with B along the z direction, 

/,(i) = Bej^o ^ /, /.(i) = /y(i) = (13) 



In this case, the ground state of Eq.([TTD will consist of a periodic structure of straight vortex 
lines, aligned with B in the z direction, with density /. 

The interaction potential G in Eq.(|TT|), is the lattice Green's function which solves 

D,p{r; - fk) = -5i,k (14) 

where is the lattice Laplacian. G{r) is most easily expressed in terms of its Fourier transform 

k 

where is the number of sites in the mesh, and 

Gu = l/K^ = l/{Kl + Kl), 

= 2 — 2 cos k^a, K'j_ = A — 2 cos k^-a — 2 cos kyU (16) 

Note K"^ ~ (/ca)^ for small ka, however they are different at large ka due to the discreteness of 
the mesh. For large r ^ a, G{f) ~ l/(47rr). 

We thus see that the interaction between vortex lines in our model is the Coulomb interaction, 



Gk = 1/K . For comparison, the mean field London interaction in the continuum[^, obtained 



by minimizing the free energy Eq.([l|) with respect to a spatially varying b{r), is Q^°^<^^ = i^(^f^'^-^ 
A^^), which is a screened Coulomb interaction with screening length A. Thus our approximation 
of taking a uniform magnetic induction b{r) = B, and ignoring spatial variations and fluctuations 
in b, corresponds formally to making this screening length infinite. Our model therefore strongly 
suppresses long length scale (small k) vortex line density fluctuations, as compared to the true 
superconductor. However we expect that the energetics of fluctuations on small length scales 
k~^ ^ A are well approximated. Provided <^ A, the model should therefore correctly describe 
the high k ~ l/a^, fluctuations describing collisions between near neighbor vortex lines, which 
are believed responsible for the melting of the vortex line lattice 0. 

It will be convenient for later to write Eq.([TTD in terms of Fourier transforms. Defining 

n,{k)^Y.M^^"'''' (17) 

i 

we have as the Hamiltonian 

FAnM = ^ E n,{k)n,{-k)G, (18) 

The k = contribution to would be [n^(A; = 0) — N f^]'^Gk=o, where is as in Eg. (PI^D . 
However as Gk=o diverges (see Eq.(p!6|)), the requirement that the total free energy be finite, 
insures that we only have configurations where the average vortex density in direction /i is 
hence the k = piece vanishes. 
III. Monte Carlo Simulations 

To investigate the behavior of the model, we carry out Metropolis Monte Carlo simulations, 
of the Hamiltonian (||), on a cubic mesh with N = L^L^ sites and periodic boundary conditions 
in all directions. Our simulations are carried out directly in terms of the phase variables 6i, with 
no additional presumptions concerning the behavior of vortex lines, beyond those implicit in 
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Eq.(^). Simulations of this same model, only working directing in terms of the vortex variables 



and using the Hamiltonian (11), have recently been reported by Cavalcanti et a/p3| 



Previously we have reported [0, |T3| on our simulations using a density of field induced vortex 
lines of / = 1/5. In this work we report on simulations using a density of lines / = 1/25, with 
an average spacing between lines of ay = 5a. This value of / was chosen so as to increase the 
ratio a^/a, and hence reduce effects due to the discreteness of our numerical mesh, while making 
/ large enough so as to have a large number of lines contained within a system of computational 
size. Our present simulations are carried out for a mesh of size L± = 50 and = 24, giving 
100 vortex lines in the ground state. The ground state of our system, we believe to be the 
periodic extension of the unit cell of straight vortex lines, shown in Fig.l. Although we have no 
way to prove that this is the true ground state, anealing with Monte Carlo simulations failed 
to result in any state with lower energy. That our presumed ground state is a square lattice of 
lines, instead of the triangular lattice expected in a continuum model, is due to the competition 
between vortex line interactions and the cubic discreteness of our numerical mesh. Most of our 
calculations represent runs of 2, 000 Monte Carlo sweeps to equilibrate, with 20, 000 sweeps for 
computing averages. Each Monte Carlo sweep refers to one pass through the entire cubic mesh. 
Errors are estimated either by block averaging, or from comparing independent runs. Henceforth 
we measure length in units where C,o = o- = ^, and temperature in units where fc^ = 1. 

A. Helicity Modulus 

The most convenient method to detect phase coherence, ie. superconductivity, in the model 



Hamiltonian (0), is to compute the helicity modulus [El]. The helicity modulus T„ in the direction 



fi, is defined as follows] 25]. Imagine applying boundary conditions such that the phase across the 
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system in the direction fi is twisted by an amount L^5. If fi = x for example, we would require, 

9{x = L^,y,z) - 9{x = 0,y,z) = LJ (19) 

Now transform to the new set of variables 

6', = 6,- in ■ (20) 

If the 9i obey the twisted boundary conditions ([19D , then these new variables 9[ are easily seen 
to obey periodic boundary conditions. In terms of the the system with the twisted boundary 
conditions is described by the Hamiltonian 



8] = JoY. ^(^: - - - (% ■ m (21) 

where ej,- is the unit vector from site i to site j. One can now form the partition function 



Z(5) = ^e-^[^-^l/^ (22) 

where the sum is only over configurations {9[} satisfying periodic boundary conditions. The 
total free energy is then 

jr(5) = -TlnZ((5) (23) 
and the helicity modulus is now defined as 

1 d'^T 



(24) 



where = LzL\. measures the stiffness to twisting the phase up in the fi direction, and may 
be thought of as the renormalized coupling Jg on infinitely long length scales. Equivalently, T 
may be viewed as the effective superfluid density, or via Eq.(|^), as determining the renormalized 
magnetic penetration length, T ~ V-^l?- When > there is phase coherence in the fi 
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direction, and when = 0, phase coherence is lost. The vanishing of in a particular 
direction therefore signals a superconducting to normal transition. 

In terms of the definitions Eg. (pT|— , we can write the helicity modulus as the correlation 
function, 

= ^ (e ^"(^^ - - ■ f^y) (25) 



•^0 



- 2, 



TN 

which can be evaluated directly in the Monte Carlo simulation. V and V" are the first and 
second derivatives of the Villain function defined in Eq.(|TO|). 

In Fig. 2 we show our results for and Tj_ = \{'^x + ^y)i as a function of temperature. 
Both heating from the ground state (Fig.l), and cooling from random are shown. The results 
clearly indicate two sharp transitions. As T is increased, there is first a transition Tci_ — 1.3Jo at 
which — * but > 0. Phase coherence is lost in the directions perpendicular to the applied 
magnetic field, however phase coherence persists in the direction parallel to the magnetic field. 
Then, at a higher T^z — 2.5Jq, and phase coherence is lost in the parallel direction as 

well, and we have the completely normal phase. In the intermediate phase, Tc± < T < T^z, the 
linear resistivity pz parallel to B should vanish, however the linear resistivity p± perpendicular 
to B should remain finite. Such an intermediate phase is in complete agreement with recent 



predictions by Feig'lman and co-workers||26||, who used arguments based on the "2(i boson" 
approximation. At the lower transition, Tc±, a slight hysteresis is seen, comparing heating versus 
cooling. No hysteresis is seen at T^z- 
B. Vortex Line Fluctuations 

We now seek to describe the nature of vortex line fluctuations in the three distinct phases 
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discussed above. We will find that the high temperature phase is one in which vortex lines wander, 
entangle, and cut through each other easily. In this phase, we also find that the density of vortex 
lines is greatly increased by the excitation of many closed vortex loops. The intermediate phase 
is one in which vortex lines are unentangled, remaining on average straight and aligned with B, 
however they are in a liquid state as far as their correlations in planes perpendicular to B are 
concerned. 

As our simulations are carried out in terms of the variables 6i, we first must locate the 
positions of the vortex lines as they thread their way through the system. We use the following 
algorithm. For each plaquette a in the system, we compute 

2vr/a + E l^i - - = 27rna (26) 

a 

where the sum is taken counterclockwise around the bounds of the plaquette, fa is the magnetic 
fiux through the plaquette as given by Eqs.(|12|,|13D, and the [Oi — Oj — Aij] are measured restricted 



to the interval [— vr, tt]. is then the integer vorticity in the 6i going around this plaquette. As 
vorticity is conserved, each unit cell of the cubic mesh which has one plaquette with a vortex 
entering the cell, has another plaquette with a vortex leaving. By computing the vorticity 
entering and leaving the unit cells of the mesh, we are able to trace out the paths of the vortex 
lines. In the event that two vortex lines enter and leave the same unit cell, we randomly assign 
which exiting segment is connected to which entering segment, for the purpose of identifying the 
paths of those lines. 

(i) Upper Transition Tcz 

One property that has received much theoretical attention, has been the entanglement of 
vortex lines, once the ground state vortex line lattice has melted. Interesting dynamical effects 
have been suggested assuming a large energy barrier exists against vortex lines cutting through 
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each other m, ^, In such a case, entanglement would provide a topological constraint 



limiting vortex line motion. A single pinned line would pin the many lines with which it is 
entangled, substantially reducing the flux flow resistance. 

In our model, it is simple to check for entanglement. We follow a prescription proposed 
by Nelson[^ 0. Let {r±i{0),f±j{0), ...} be the positions of the vortex lines in the xy plane, 
as they pierce the plane at z = 0. Let {r±i{Lz),r_Lj{Lz), ...} be the corresponding locations of 
these vortex lines as they pierce the z = plane. Since our simulation uses periodic boundary 
conditions in the z direction, parallel to B, all the f±i{Lz) must connect onto some r±j{0), so 
that the vortex lines remain continuous in the periodically extended system. If r±i{Lz) connects 
to r±i{0), we say that the line i reconnects to itself. If we view the periodic boundary condition 
along the z direction as representing the circumference of a three dimensional torus, then line i 
winds only once around the torus, before reconnecting to itself. In contrast, if f±i{Lz) connects 
to r±j{0), j 7^ i, then lines i and j are part of a larger multiple connection of lines that winds 
more than once around torus, before reconnecting to itself (see Fig. 3). For an entangled vortex 
line state, we would expect there to be many such multiple connections, while for disentangled 
approximately straight lines, there should be only simple self-reconnections. In Fig. 2, we plot 
the fraction of lines which reconnect to themselves, R, versus temperature, which we compute 
while cooling the system from an initial random configuration. R drops sharply to zero at T),^, 
from its low temperature value of one. We thus identify this upper transition, where phase 

— * 

coherence in the direction parallel to B is lost, as a transition from straight disentangled lines, to 
wandering entangled lines. Further, the fact that the multiply connected entangled lines present 
at high temperature, completely disentanlge to simple self-reconnected lines below Tcz, suggests 
that the lines may easily cut through each other as they disentangle. 
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As an alternative criterion for vortex line entanlgement, we can consider the transverse fluc- 
tuations of the vortex lines {v?{lz)) as a function of length < < parallel to 5, 

{n\Q) = -^j:T.{[rU^)-rUzo)r), Zo = IJ2 (27) 

where f\Lj(-2) is the position of line i in the xy plane at z, Ny = fL\ is the number of magnetic 
field induced lines, and the sum is along the lengths of these lines for a distance Iz- Vortex 
lines in our model may in general consist of two types: the magnetic field induced lines which 
thread the entire system from z = to z = L^, and thermally excited closed vortex loops. In 
the evaluation of Eq. ( PTj ) we do not include sums over vortex line segments which belong to the 
thermal loops. In the event that a loop intersects a field induced line, we randomly assign which 
line segment "belongs" to the field induced line, and which segment "belongs" to the thermally 
excited loop. 

For straight, unentangled vortex lines, we expect that (u^) will saturate to a constant as 1^ 
gets large. For wandering, entangled vortex lines we expect [|l| that (m^) will grow with increasing 

Iz as, 

{u\h)) = ul{T) + D{T)h (28) 

D{T) measures the diffusion of the line in the transverse direction, as one moves down the z 
axis, and it can be used to determine the "entanglement correlation length," |I| 

UT) ^ (29) 

measures the distance one has to move in the z direction, in order to have a transverse 
fluctuation equal to half the average spacing between the vortex lines. is the length scale 
on which the fluctuating lines can approach each other, leading to a mutual winding and hence 
entangling. 
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In Figs. 4a, 6 we plot {u'^{lz)), versus Iz, for several different values of temperature. For T > 
2.5 Jq ~ Tcz (Fig.4a), the data points for the three largest values of 1^ show a clear linear increase 
with Iz, as expected from Eq.(]28|) for an entangled phase. For T < 2.5Jo (Fig.46), there remains a 
systematic curvature in the data, suggesting a saturation of (u^) at large Iz, as in an unentangled 
phase; the data here is limited by finite size effects, and runs in a system with larger Lz would 
be needed in order to clarify the situation. Obtaining D{T) from a fit of Eq. (|28D to the three 
largest values of lz, for the data of Fig.4, we compute the resulting C,z from Eq.(^) and show the 
result in Fig. 5. 

In a similar manner, we can compute the "cutting correlation length," C,c, which we define 
as the average distance one has to move in the z direction, between two successive cuttings of a 
given vortex line. Within our model, we define a cut as whenever there are two or more vortex 
lines entering and leaving a given unit cell of the numerical mesh. Such a situation represents 
an overlapping of the vortex cores, and there is no constraint whatever on how the lines will 
subsequently fluctuate out of this unit cell; lines are free to cut through each other, or even 
to cut and reconnect different segments of different lines. If Nc is the average number of total 
cuttings observed in the simulation, then 

UT) ^ ^ (30) 

We plot i^{T) in Fig.5. 

Comparing the entanglement and cutting correlation lengths, and ^c, we see that for 
T > 2.5Jo — Tcz, C,z — ^c- Thus, entanglement and cutting are occurring on the same length 
scales for the region above T^z, and here we would expect entanglement to provide little constraint 
on the motion of vortex lines. Below Tcz, we see growing faster than ^z- However, as discussed 
above in connection with Fig.46, we believe that our data here is limited by finite size effects, and 
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hence the values shown in Fig. 5 are not a rehab le indication of behavior in an infinite sample. It 
therefore remains unclear whether there exists some narrow range of temperature, at or below 
Tcz, in which entanglement might lead to constraints on vortex line motion. 

Finally, as an additional measure of the fiuctuations of the vortex lines, we calculate the 

— * 

average density of vortex lines cutting the planes perpendicular to B, 

n^^l^Ein'M (31) 

i 

where the sum is over all dual sites i, and nz{i) is the vorticity passing through the plaquette 
in the xy plane above site i. In virtually all cases, = 0, ±1, and so Uy counts the density 
of lines, irrespective of the sign of the vorticity. If the only vortex lines present in the system 
are the magnetic field induced lines, and these lines wander through the system in a directed 
fashion (ie. all lines pierce each z plane only once; there is no "snaking back") then one has 
Uy = f. If however, there are additional closed vortex line loops excited, or lines are not directed, 
then riy will be larger than /. In Fig. 6 we plot ny{T). At low temperatures T ^ 2.0Jo, we see 
that rij, ~ /, and hence there are only directed field induced lines. However for T ^ 2.0Jo, 
grows dramatically, becoming n„ ~ 2/ at the transition Tcz- Thus, by the time the transition 
is reached, there are many thermally excited closed vortex loops in the system. This suggests 
that the transition at T^z is mediated, not by the wandering of the field induced vortex lines, 
but rather by the proliferation of thermally excited vortex line loops. The behavior of transverse 
fiuctuations and cuttings, shown in Figs. 4 and 5, may be understood as arising from the growth 
of loop excitations, that intersect and link up neighboring field induced lines, rather than the 
direct wandering of the field induced lines. In fact. Fig. 6 shows a qualitative resemblance to the 



density of vortices as one crosses the transition, in the ordinary three dimensional XY model 
(corresponding to 5=0). Similarly, a plot of the specific heat C, which we compute by numerical 
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differentiation of the average energy and show in Fig. 7, shows the same cusp-hke feature at Tcz 
as in the 3d XY model[^. We note that these loop excitations, which we find important for 
the transition at Tcz, are exphcitly excluded in the "2(i boson" and effective elastic medium 



approximations for vortex line fluctuations, as well as in other earlier simulations . 
(ii) Lower Transition Tc± 

The low temperature phase T < Tc± is characterized by an ordered vortex line lattice, with 
periodicity as given in Fig.l. The vanishing of T_|_ above Tc±, together with the result that 
~ 1 for Tc± < T < Tcz (see Fig. 2), indicates that at t Tc± the vortex line lattice melts into 
a liquid of unentangled vortex lines. One way to measure this transition, is by considering the 
local orientational order of vortices in the planes perpendicular to B. These vortices are just the 
intersections of the vortex lines, with a particular plane at constant z. We define the n-fold local 
orientational order, in terms of the quantity 

1 1 

*n = 7]vE7zEe^"'^^^^"^^^^ (32) 

z,i * j=l 

In the above, z labels the planes perpendicular to B, i labels the vortices in the plane z, and 
j labels the Mi nearest neighbor vortices of i in plane z. is the angle between the line 

from i to j, and the line from i to j + 1, which is the next neighbor after j which one encounters 
in rotating around i (see Fig. 8). The nearest neighbors of i are determined by the method of 



Voronoi polygons ||32||. For perfect square orientational order, with all = tt/2, such as 

in the ground state of Fig.l, we have \l/4 = 1, while \l/6 = e*'^. For perfect hexatic orientational 
order, with all = tt /3, = 1, while \l/4 = e*^'^/^. In Figs. 9a, 6, we plot respectively the 

amplitude, and phase of ^4 and \l/6- Both amplitudes show a clear kink at T ~ I.SJq — Tc±. At 
this temperature, the phases indicate a clear cross over from the 4-fold order of the ground state 
vortex line lattice, to the local 6-fold order of a vortex line liquid. 



For a more global picture of the transition from vortex line lattice to vortex line liquid, we 
have computed the correlation of vortices in the planes perpendicular to B. If we define the 
vortex structure function |jl], |^ as 

S{k^,h) = ^Mk)n,i-k)) (33) 

with k = {kj_,kz) and nz{k) as in Eq.(|r^), then the correlation between vortices in different 
planes separated by a height z, is given by the Fourier transform of Eq.(|53D with respect to kz, 

S{k^,z) = ^J2Sik^,h)e"'^' (34) 

The correlation between vortices in the same plane is then just S{k±,z = 0). In Figs.lOa-d 
we show intensity plots of S{k±,z = 0) for k^^ky in the interval [— 7r,7r], for several values of 
temperature, upon heating the system. At T = 0.6Jo < Tc±_ (Fig. 10a) we see sharp Bragg peaks, 
corresponding to the periodic vortex line lattice of Fig.l, superimposed on a small continuous 
background. The width of these Bragg peaks, corresponds to the finite resolution of wavevectors 
allowed by our finite system, A/c^ = In jL^. At T = 1.5 ^ Tc^ (Fig. 106), we see an approximately 
rotationally invariant structure, with peaks of finite width. The first peak is at a radius fc_L = 
lixja^^ corresponding to the average separation between vortices; the smaller, second peak, is 
at /c_L = Air/ay. Such a structure function is characteristic of a liquid. The lack of complete 
rotational symmetry in the peaks (aside from a slight 4-fold distortion to be expected from 
the underlying cubic mesh) we believe is due to the finite time duration of our simulation; 
we see no evidence for long range orientational order as might be present in a hexatic line 
liquid[^]. As T increases (Fig. 10c), the second peak disappears into the continuous background. 
At T = 2. 5 Jo — Tcz (Fig.lOc?), both peaks have virtually disappeared, and only a continuous 
background remains. Thus for T ^ Tcz-, even the short range order, characterized by the average 
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spacing between field induced vortex lines, is lost; this is consistent with the proliferation of 
many vortex line loops, as evidenced by Fig.6. 

In the intermediate phase, we find that the system relaxes very slowly to equilibrium. We 
illustrate this by computing the time dependent correlation, 

9{t) = ^E{K(MK(i,0)) - (n,(i,t))(n,(i,0))} (35) 

i 

where the time t is just the number of Monte Carlo sweeps through the entire cubic mesh. We 
plot g{t) at T = 1.9Jo in Fig. 11. After an initial fast decay, we observe a slow algebraic decay, 
g{t) ~ t^'~- with ~ 0.23. We believe that it is this very slow decay which is responsible for the 
small but finite values of that are seen for T > T^^ in Fig. 2, rather than usual finite size 
effects. When we doubled the number of Monte Carlo steps in our simulation, we observed 
in this region to decrease. 

Because of this slow equilibration, as the system is cooled back below Tc_l, we have been unable 
to regain the ordered vortex line lattice of Fig.l. Instead, we cool into a structure consisting of 
perfectly straight lines, ordered into local domains of the ground state, with frozen in domain 
walls. In Fig. 12 we show an intensity plot of S{k^,z = 0) at T = O.GJq, for such a cooled 
configuration. The two rings of peak intensity which exist for T > T(.± (see Fig. 106), are now 
narrow and "squared off." This feature of freezing in domain walls, as well as the hysteresis in 



seen in Fig. 2, suggest that the transition at Tc_l may be first order Another possibility is that 
is not a true thermodynamic transition at all, but only a cross-over depinning temperature 
where the effective energy barriers to vortex line motion grow so large, that vortex motion is 
frozen out on the scales of our simulation. 
C. Vortex Structure Function 

In Fig. 10, we used the structure function S{k±,kz) to qualitatively determine the nature of 
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the vortex line structure in the different thermodynamic phases. Here we present a more detailed 
analysis of S{k±, k^), from which we can extract information about the effective elastic constants 
of the model. 



Following the works of Marchetti|]35| and Nelson and LeDoussal0, the natural form for the 
coarse-grain averaged free energy of our model, valid on hydrodynamic length scales fcj^ > a„ 
in the vortex line liquid phase, is 



'^^^ J k \ ^J.±=x,y 



n^^{-k)\+Q{k)6n,{k)6n,{k)j (36) 

where the 6n^{k) give the deviation from the average uniform density and are now viewed 
as independent continuous variables. C44 and c/ are the tilt and bulk moduli of the vortex lines 
on the numerical mesh. 

Using the constraint that vortex lines must be continuous, ie. the discrete divergence of ^^(i) 
vanishes, 

^K(i + /i)-n^(i)] = 0, or ^n^(^)[e^^.-l]=0 (37) 

M k 

one can do the appropriate Gaussian integrals^ over 6n^{k) and compute the structure function 
of Eq.( 



S(k^, k,) = ^ ^ (38) 

Q(fc)Ki + C44(fc)K2 

where is related to A;^ as in Eq.(|T6|). The factor Kj_ in the numerator of the expression above, 
determines the qualitative form of the continuous background of S{k±, z = 0) seen in Figs. 10. 
From Eq.(|36|) and the form of the bare vortex line interaction given in Eqs.(|l^,0), we expect 

ci{k) = 47rVViGfc, C44(fc) = ci{k) + fei{k,) (39) 

where we write Ji instead of Jo to allow for possible renormalization of this coupling in the 
coarse-graining procedure, ei is the single vortex line tension, which depends on kz, due to the 
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interaction between the different segments of the single hne. The difference C44 — q just gives the 
additional energy needed to create the elongation of the vortex lines described by the transverse 



components of the fluctuation Substituting Eg. ([39|) into Eq.(|38|) gives 



TKl/Sik^, h) = Att'J, + [ei{h)/f]K^, (40) 

We now compute S{k±,kz) directly within the Monte Carlo simulations. Calculations are 
done for values of kz = 27rm^/L2, rriz = 0, 1, and for the special values of k± = k±x 
and kj_y, where kj_ = 27rm_|_/L_|_, m_|_ = 1, L_|_/2. By symmetry, we expect that S{kj_x,kz) = 
S{k±y, kz), hence we combine our results by defining S{k±, kz) as the average of S in these two 
special directions. 

In Figs. 13a — /, we plot TKj_/S, as a function of Kj_ and K"^ for temperatures T = 1.0 Jq 
(in the vortex line lattice phase), T = 2.25Jo (in the intermediate phase), and T = 2.75Jo (in 
the normal phase). In the line lattice phase (Figs. 13a, 6), we see that TKj_/S is approximately 
independent of k±, obeying the hydrodynamic form of Eq.(^OD at all k. In the intermediate, 
unentangled line liquid phase (Figs. 13c, d), a strong k± dependence appears, with a minimum 
appearing at Kj_ ~ 1.2, very close to the value of 1.38 which corresponds to k± = 271 /ay. In the 
normal phase (Figs.l3e, /) as T increases, this minimum decreases to smaller k±, reflecting the 
increasing density of vortex lines (see Fig. 6; this also explains why the minimum in Fig. 13c is 
not exactly given by k± = 2'7r/a^). This minimum as a function of kj_, is a reflection of the peak 
in S{k±,kz) which occurs isotropically in the k± plane at ^ 27r/a^, in the line liquid phase. 
It may be viewed as an effective softening of the bulk modulus C44(/c) at these large values of k±. 
In the line lattice phase, this isotropic peak contracts into discrete Bragg peaks, none of which 
lie along the x ot y directions (see Fig. 10a). Hence the kj_ dependence of TK\/S flattens out. 

Although Eqs. (|38|-^) were proposed for the line liquid phase, similar results have been given 
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for the line lattice phase, using an effective elastic approximation |p. In this case the only change 
is that one must replace q in Eq.(p8D with cn = Q + Cqq, where Cqq is the shear modulus. This 
results in the addition of the term cee-ft'l to Eq.(pOD. It is unclear whether or not the differences 
observed in the k±^ dependence of Figs. 13a and c, may be interpreted in terms of the vanishing 



of the shear modulus ||33| at the melting transition Tc±. 

Comparing the curves of TK\/S versus Kl, for different fc^, in Figs.l3(i and /, one sees 
that the dependence of TKj_/S on k± and fc^ cannot be separated out into the simple form 
4:71^ Ji{k^) + [ei{kz) / f]K'^ as might be expected from the "2(i boson" approximation 0. As we 



have noted in earlier work||13||, the more complex dependence we observe, may be interpreted as 
adding a k± dependence to the line tension ei. Using Eq.(^) to define the effective line tension 
as k± — > 0, 

/ 



eiik,) - ^2 



TKi TKj 



(41) 



5(0, A;,) 5(0, 0)_ 

we plot in Fig. 14 ei{kz) versus Ini^^^ several temperatures T. The plots approach straight lines 
at large K^. This is in qualitative agreement with the expression for the line tension obtained 
from an elastic approximation, expanding about the line lattice in a layered superconductor 0], 



2 



eJk,) = , ° In 
2(47rA)2 



2 



.,max 



(42) 



IP- -I- K'^ 

where for the line lattice, k± fnax ~ ^r/^O) a-^d for A ^ a^, k^^min ~ vr/a„. For our data, we find 
a reasonable fit with k±^max ~ (in our model = 1), but k±^rnin ~ 0. 

From the above analysis, we can extract several parameters of interest. Considering the {k±, 
kz) — *■ limits, we fit our data to the form of Eq. (|40|) and extract the values of the hydrodynamic 
parameters Ji and ei(0), which we plot versus T in Figs. 15a and b, respectively. Defining the 
prefactor of the logarithm in Eq.(^2[) as (7r/2)J2 (as motivated by Eq.(^), we can fit the large 
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part of the data in Fig. 14 to ~ hi{Kl) and plot the resulting J2 in Fig. 15a. We see that 
Ji and J2 are of comparable magnitude, except at the very lowest temperature T = O.GJq, and 
above Tcz — 2.5Jo, where J2 decays to zero, while Ji remains approximately constant. 

Finally, we can extract a correlation length along the field direction, by computing the Fourier 
transform of S{k±,kz) with respect to k^, as in Eg. (|3^ . Fitting the resulting S{k±,z), at large 
z, to the form 

g-^/?{fcx) _^ g-(i.-2)/^(fe±)l (^43) 

(this form being chosen so as to account for the periodic boundary conditions along z) , we plot the 
correlation length C,z = ^{k± = 27r/at,) versus T in Fig. 5. For T < 2. 25 Jo we found no measurable 
decay in S{k±, z). For T > 2.25 Jo we find ~ the cutting correlation length. This reinforces 
our earlier observation that cutting and entangling appear to occur on the same length scales. 
We can use the value of ^z just obtained, to estimate the hydrodynamic line tension. For finite 
^zi we expect that the effective fc^ ^ line tension can be obtained from Eq. (|4^) (appropriate for 
the line lattice where ^z 06) by imposing a small kz cutoff at fcz ~ vr/.^^. Hence we estimate, 

ei (0) ~ ei ikz - Tx/lz) ~ vr J2 In (44) 

We plot the result in Fig. 156, where we have set = Lz for all T < T^z- Although we observe 
quite reasonable agreement with the hydrodynamic result obtained from a direct fit to Eq.(^Dp, 
we note that the dominant temperature dependence comes from the behavior of J2. We do not 
have large enough of a range of ^z to provide a sensitive test of the In^^ dependence. 
IV. Conclusions and Discussion 

To conclude, we have studied a model of fluctuating vortex lines, and find three distinct 
thermodynamic phases. Upon heating, the system has sharp transitions from a vortex line 
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lattice, to an unentangled vortex line liquid, to an entangled vortex line liquid with many closed 
vortex line loops and much vortex line cutting. The lower transition is associated with the loss 
of phase coherence in the superconducting order parameter, in directions perpendicular to the 
applied magnetic field. The upper transition is associated with the loss of phase coherence in the 



parallel direction. This behavior is consistent with that argued by Fiegel'man and coworkers p6 
Similar behavior has been argued by Glazman and Koshelevf^, for low magnetic fields B < 
Bcr^^- Nelson[0] originally suggested an unentangled fine hquid phase as a finite size effect. 



when C,z ~ Lz- We cannot rule out the possibility that the difference we observe between Tc± 
and Tcz is such a finite size effect. However, if this is so, our observation that Tcz — ^Tc^ 
would indicate a rather large critical region about the single transition of the infinite system. 
Similarly, more detailed finite size studies would be needed for one to clarify the nature of the 



phase transitions in our model. In our earlier work|T^, |T3[, at the higher density / = 1/5, 
we found evidence for only a single transition. It may be that there is indeed only a single 
transition at higher densities, or it may be that the width of the intermediate region becomes 
too narrow to distinguish. It is interesting to note that in a related model, the 3dXY model 
of stacked triangular planes with antiferromagnetic coupling within the triangular planes and 
ferromagnetic coupling between the planes, which has the same symmetry [p^j as our model with 



/ = 1/2, there is only a single transition, however it is a tetracritical point p8| 



Our model has involved two key approximations in its mapping to a realistic superconductor. 
One is the discretization of space into a cubic mesh, on the length scale of the bare coherence 
length ^0- The discretization along the z axis may be viewed as representing the layered structure 



of the high Tc materials, as in a Lawrence- Doniach model||3^. The discretization in the xy plane 



however, represents an artificial periodic pinning potential for vortex line fiuctuations. Our hope 
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was that this pinning potential primarily acts to remove the mode associated with a uniform 
translation of the entire vortex line lattice. However the "forces" associated with this pinning 
potential are strong enough to tilt the balance in favor of a square line lattice in the ground 
state, as opposed to the triangular lattice expected in a continuum. One may worry that these 
forces also significantly suppress vortex line fluctuations, leading to a stabilization of the line 
lattice, as well as phase coherence of the superconducting order parameter, at low temperatures. 
However, in Figs. 15 we see a sharp drop in Ji and ei(0) at a T < I.OJq) noticeably below the 
lower transition Tc± ~ I.SJq- If we interpret T as the energy scale associated with the energy 
barriers of the periodic pinning potential, then one expects that the transition at Tc± is a true 
melting transition of the line lattice, rather than just a depinning of the lines. 

The second approximation has been the assumption of a uniform magnetic induction B inside 
the superconductor, ie. the A ^ oo limit. We first discuss the effects of this approximation on 
vortex line fluctuations. As mentioned in Section II, setting A — > oo suppresses long wavelength 
(small k < A^^) density fluctuations, while being a good approximation for high k ^ A^^ 
fluctuations. It is generally believed0, |^ however, that it is high k± ^ 27r/a„ shear fluctuations, 
which do not alter the long wavelength density, that are responsible for melting of the vortex 
line lattice. Ikeda et a/.[^ specifically note that the melting transition, as determined from a 



Lindemann criterion, is largely unaffected when comparing the finite A to the A ^ oo cases, 
providing A a„. Hence we expect behavior in the line lattice phase, and the melting at T^x, to 

— * 

be largely independent of this uniform B approximation. Similarly, in the normal, entangled line 
liquid phase, T > T^z, one expects this approximation to be good for describing local properties 
(such as entangling, cutting, and S{k) for k > A^^) provided A ^ ,^2, the correlation length along 
the z axis. In this case, fluctuations impose a shorter cutoff on effective interactions than does 
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magnetic screening. It is less clear how valid the approximation is in the intermediate phase, 
where .^^ — oo, and there is no other natural length scale along the z direction with which to 
compare A. However we note that Feigel'man|26[| predicts a similar unentangled phase, working 



within the "2rf boson" approximation, in which the interaction between vortex lines along the z 
direction is short ranged. Further work will be necessary to determine whether the unentangled 
line liquid phase observed here, is stable against fluctuations of B. 

The effect of the A ^ oo approximation on fluctuations of the phase of the superconducting 
order parameter, and the existence of true long range order in the gauge invariant order parameter 
correlation function, has been a subject of much discussion^, 1^, 0, 0. Moore|^ and Ikeda 



et a/. [40| have shown that in three dimensions, the gauge invariant order parameter correlations 



decay algebraically, in the line lattice phase, for the \ ^ oo case. However for the finite A case. 



they decay exponentially. Such algebraic decay is consistent with a finite helicity modulus 
T, as we observe in our model. Exponential decay, for finite A, would imply T ^ 0. However as 
shown in Refs. the length scale for such exponential decay is given by (7rJo/T)AK, and 

is hence comparable to the size of macroscopic samples, for the large k high Tc materials. Thus 



the A — s> oo limit should be appropriate in such situations p4|. 

Although we have not attempted any direct mapping between our model and a specific high 
Tc material (we took isotropic couplings, and ignored all T dependence of our bare parameter Jq), 
it is tempting to try to discuss resistivity measurements in very pure YBCO single crystals [|7|, |i5[|, 
in terms of the behavior we observe in our model. Above the upper transition T^z-, where there 
are many vortex line loops, much wandering and entangling of the field induced vortex lines, 
and easy cutting of vortex lines, we might expect substantial flux flow resistance due to the 
flow of effectively independent vortex line segments. This region possibly describes much of the 
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broadening of the resistive transition, in comparison to the B = case, and could be explained 



with traditional flux flow models]^. Resistivity, p^, for current / || B, should be comparable 
to p± for I J- B (apart from effects due to explicit anisotropy of couplings), as there are almost 
as many vortex line segments oriented perpendicular to B as parallel to B (from closed loops, 
and transverse segments of wandering lines). As T is cooled towards Tcz we would expect 
Pz to decrease substantially compared to p±, as the vortex lines disentangle, and straighten 
parallel to B, and there are fewer loop excitations. Such behavior is consistent with that seen 



in experiments! 45 1 where B is applied in the ab (CuO) planes, and hence there is no explicit 
anisotropy effects for p^ versus p±. Below T^z, the linear resistivity pz should vanish, however 
the linear part of p± would remain finite. The transition at T^z might offer an explanation for 
the "kink" observed in p± in the tail of the resistance curve, for very pure samples with B along 
the c axis]^. Finally, p± would vanish below Tc±. 
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Figure Captions 

Figure 1. Ground state vortex line lattice for a magnetic induction oi f — B^^/^q — 
1/25 flux quantum per unit cell of the numerical mesh. The view is along the direction 
of B, and (+) locates the positions of the straight vortex hues. 

— * — * 

Figure 2. Hehcity modulus along direction of B, and T_l perpendicular to B, 
for both heating and cooling. The vanishing of T;^^^ indicate two separate phase 
transitions. R is the fraction of vortex lines which reconnect to themselves, when the 
system is periodically extended along the z direction (see section IIIB(i)). Vanishing 
of R indicates entanglement of lines. 

Figure 3. Schematic, showing how vortex lines must connect onto each other, due 
to the periodic boundary conditions along the z axis. Solid lines represent vortex 
hues. Dashed lines indicate how the ends of the vortex hues must match up on the 
z — and z — planes. 

Figure 4. Transverse fluctuation of vortex lines {u^) as a function of length along the 
line for several different temperatures. Linear growth with 1^ indicates wandering, 
entangled lines. Saturation with increasing 1^ indicates unentangled lines. 

Figure 5. Correlation lengths along the direction of is the "entanglement" 

length, determined by the transverse fluctuations of Fig. 4. is the "cutting" length, 
giving the average distance between two cuttings (intersections) along a single vortex 
line, is the decay correlation length as extracted from the vortex line structure 

— * 

function, S{k_i_,z) (see section IIIC). 
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Figure 6. Average density of vortex lines n^, piercing planes of constant z. At low 
T, is the magnetic field induced density of / = 1/25. increases with T due to 
the excitation of closed vortex line loops. 

Figure 7. Specific heat C versus T. 

Figure 8. Schematic, showing the measurement of local bond angles used 
for the calculation of local orientational order in Eg. (^2]). 

Figure 9. 4-fold and 6-fold local orientational order parameters, ^4 = |\l/4|e*''* and 
\l/6 = |^&4|e'^''. The kink in |\&4,6| indicates Tc± ~ 1.3Jo, while the crossover in -^^^ 
indicates a transition from local 4-fold to 6-fold order. 

Figure 10. Intensity plots of vortex structure function S{k±, z = 0), for various tem- 
peratures, kx is the horizontal, and ky is the vertical direction. Brightness measures 
the magnitude of S, as in a diffraction pattern. The 5-function peak at fc^ = is the 
temperature independent constant 5(0, z = 0) = pL\, which sets the relative inten- 
sity scale, (a) shows a lattice of Bragg peaks, characteristic of a vortex line lattice, at 
T < Tc±. [b) and (c) are characteristic of a vortex line liquid, for Tc± < T < T^z- (d) 
is for T > Tcz, and the dominating continuous background indicates the proliferation 
of closed vortex line loop excitations. 

Figure 11. Time decay of real space vortex correlation function of Eq.(|35|) for 
T = 1.9 Jo, in the intermediate phase. The linear behavior on the log-log plot indicates 
algebraic decay with the power 0.23. The unit of time is one Monte Carlo sweep 
through the entire numerical mesh. The solid line is a guide to the eye. 
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Figure 12. Intensity plot of vortex structure function S{k±,z = 0), for T = O.6J0, 
upon cooling the system through Tc±. 
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Figure 13. Structure function data, plotted as TKj_/ S{k±, kz), versus Kf^ and K. 
for several different temperatures. are related to k^ as in Eg. ([T^) . The different 
curves versus Kj_, are for different values of kz = 2'Kmz/Lz. The different curves 
versus K^, are for different values of k± = 27im±/L± {m± = 10 corresponds to 
k± = 27r/a^). (a) and (6) are for T = 1.0 Jq in the line lattice phase; (c) and [d) are 
for T = 2.25 Jo in the intermediate, unentangled line liquid phase; (e) and (/) are for 
T = 2.75J0 in the normal, entangled line liquid phase. 

Figure 14. Line tension ei{kz), from Eq. (PD , using data as in Fig. 13. The linear 
dependence on InK"^, at large K"^, agrees quahtatively with the theoretical prediction 
of Eq.(|2D. is related to k^ as in Eq.(|16|). 



Figure 15. (a) Coupling Ji, and [b) line tension ei{kz 0), obtained from a fit of 
TK\IS to the hydrodynamic limit, Eq.(|^). (a) Coupling J2, obtained from a high 
fit of the line tension data of Fig. 14, to the form Eg . (|4^) . (6) Hydrodynamic line 



tension 7172 In obtained from Eq . (|4^) , using structure function correlation length 
iz of Fig.5. 



36 



